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X. INTRODUCTION 


The general goal of this grant is to undertake research on „the gener- 
ation of tunable visible, infrared, and ultraviolet light, and on the con- 
trol of this light by means of novel mode-locking and modulation techniques. 
During this reporting period the following projects have been active: (l) 

transient mode-locking of the Nd:YAG laser and generation of short tunable 
pulses in the visible; (2) investigations of the alkali metal inert gas 
excimer laser systems; (3) techniques for frequency conversion of high 
power and high energy laser radiation; and ( 4 ) high average power blue and 
UV laser light sources. These projects are described in detail in Section 
II of this report. Although the coherent vacuum ultraviolet holography pro- 
ject is currently inactive, a paper based on that work has been submitted 
for publication. It is entitled "Effects of Focusing on Third Order Non- 
linear Processes in Isotropic Media" by G. C. Bjorklund and is attached as 
Appendix A. 
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XI. SUMMARY OF PROGRESS 


A, Transient Mode-Locking of the Nd:YAG Laser, and Generation of Short, 

Tunable Pulses in the Visible 
*■ 

(D. VI . Phillion, D. J. Kuizenga, and A. E. Siegman) 

For the past several report periods this, project has been concerned 

with the theoretical understanding and experimental optimization of short- 

pulse generation in NdrYAG lasers when these lasers are continuously pumped, 
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repetitively Q-switched, and simultaneously mode-locked. J During the past 
period this work has led to the completion of a Nd:YAG laser that generates 
a continuous train of single mode-locked pulses at a repetition rate of 500 
single pulses per second, each pulse having a peak power of 0.5 -1.0 MW and 
a pulsewidth of ~ 6o psec. The technique of "pre-lasing” developed under 
this contract is used to overcome the transient build-up problem for mode- 
locking of the laser, so as to achieve this pulsewidth. A single pulse is 
selected from each burst of repetitively Q-switched and mode-locked laser 
output by means of a z-cut LiNbO^ electro-optic switch. The kilovolt level, 
nanosecond rise time pulse that activates this switch is synchronized to the 
pulse burst by a fast photodiode, fast logic circuitry, and an avalanche 
transistor chain for generating the pulse, 

- - 

Ph.D. dissertation completed, now at Lawrence Livermore Laboratories. 


The resulting train of 1.06** jjm pulses can then be doubled to 532 ran 

* 

■with > 35^ conversion efficiency in an elliptically focused X,iIO^ crystal; 
and we are presently working on further doubling of the pulses to 266 ran 
in a CDA crystal. 

56. 

This laser has a peak power (lCr - 10 W) much closer to a flash-pumped 
. 7 Q 

mode-locked YAG laser (10 - Kr W) than to a cw-pumped mode-locked YAG 
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laser (10 - 10 W), at a substantially higher repetition rate than in a 

flash-pumped laser. Yet it has essentially the same stability, spectral 
purity and reproducibility as does a cw mode-locked laser. Hence it will 
be ideally suited as a high-repetition-rate source for many short-pulse 
applications, including as a master oscillator for lunar ranging experi- 
ments; for airborne optical radar profiling of water turbidity; for non- 
linear optics experiments; and for basic physical measurements of fast 
atomic and molecular relaxation processes, an area we are presently pursuing. 

The pulsewidth of this laser is, however, still substantially wider 
than the limiting value determined by the YAG atomic linewidth. Active 
mode locking as used in this laser in general does not lead to a completely 
linewidth -limited pulse. A more dynamic form of mode-locking, such as is 
provided by a saturable absorber, is required to achieve the limiting 
pulsewidth. Pulsewidths as narrow as ~ 15 psec should be possible in a 
linewidth -limited YAG laser. 

We are, therefore, now examining the possibility of adding some sort 

of passive pulse-shortening mechanism to our actively mode-locked laser. 

% 

Analytic studies have been started, to be followed later by experiments 
once theoretical predictions are completed. Two approaches appear promising. 
One ,is the straightforward use of a passive saturable absorber in combination 
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with the active modulator* This approach will undoubtedly lead to some 
further pulse shortening, but it is not yet clear whether this improvement 
will cost more than it is worth in terms of decreased peak power output. 

The second approach is the use of the optical Kerr effect to produce pulse 
self -chirping together with an intracavity dispersive element such as a 
Gires-Tournois interferometer to compress the chirped pulses. This is a 
concept that has been considered several times for single-pass pulse com- 
pression outside of high-power lasers. It does not seem to have been ex- 
plored,, however, for intracavity multi-pass pulse compression. Preliminary 
studies show that substantial chirping should be readily obtainable, and 
the pulse compression stage will be evaluated next. This approach could 
offer very significant advantages, including very small losses and a pulse 
compression that increases as the pulse narrows, rather than saturating 
with pulse growth as does a saturable absorber. 

We also intend to explore the possibilities of using the primary burst 
of pulses from our doubled laser to pump a synchronously mode-locked dye 
laser. This would give a tunable short-pulse output available in synchronism 
with the YAG pulses, a combination that could be ideal for practical applica- 
tions including, for example, airborne optical radars for mapping algae dis- 
tributions in coastal and lake waters. Substantial progress on both the 
above topics is intended during the next reporting period. 
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B , Studies Concerning the Possibility of a Li-Xe Excimer Laser 

(j. H. Newton, J. F. Young, and S. E. Harris) 

The goal of this project is to investigate the alkali metal excimer 
laser systems. In particular we have concentrated on the alkali-rare gas 
excimer s . These excimer systems may prove to be useful as high efficiency, 
high power density amplifiers or lasers in the visible or near infrared 
portion, of the spectrum. 

Excimers are molecules which exist only in the excited state. The ex- 
cited state exhibits an attractive potential while the ground state is re- 
pulsive. Physically, the mechanism for excimer laser action is as follows: 
a ground state alkali atom is excited and then combines with a ground state 
rare gas atom to form the molecule. ' When stimulated emission takes place 
the excimer drops to the ground state, dissociating into the free alkali and 
rare gas atoms. Thus the system is effectively a four-level laser, with a 
small energy defect. 

A particularly interesting property of the alkali excimer laser systems 
is that a wide variety of pumping methods may be used to excite the alkali 
atoms. Flashlamp, electric discharge, or E-beam pumping may be used to 
excite the free alkali atom. For example, a typical cross section for 
electron excitation of alkali resonance transitions is 50 which greatly 
exceeds other low-energy inelastic cross sections. Such pumping could offer 
a very efficient means of funneling energy into the upper laser level. 
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Initially we will use optical (laser) pumping which will allow us to excite 
specific levels at well known rates. 

The alkali-Xe family of exciraers are of most interest initially because 
the upper state is the most attractive relative to the other rare gases. 

This gives rise to greater population inversions and, hence, greater gain. 

Our calculations for the Na-Xe excitner indicate reasonable gain, but, un- 
fortunately, our measurements of the Na^ absorption indicate that the dimer 
absorption will dominate, preventing net gain. We have subsequently in- 
vestigated the other alkali systems and have determined that, with the ex- 
ception of the Li-Xe vs. Li^, the alkali-dimer A-X band absorption extends 
to longer wavelengths than the alkali -rare gas excimer gain. Thus, we have 
concentrated our efforts on Li-Xe. 

We have performed two different gain calculations which are in reason- 
able agreement. In the first, we calculated density distributions for the 
excited and ground states based on a Leanard-Jones fit to the excited state 
potential of Baylis^ and a uniform potential for the ground state. Popu- 
lation inversions for various pump conditions, temperatures, and pressures 
were determined. Optical cross-sections were determined from fluorescence 

data. The second calculations were based on Phelp’s method using poten- 

1 

tial curves calculated by Baylis . Predicted gain for Li-Xe is approximately 

per cm for 3*90 torx Li , 2 atm Xe, and 30$ excitation fraction. 

Figure 1 shows a schematic of a typical laser cavity. The length l 

is the gain medium and d is the dead or empty space in the cavity. The 

net gain of the laser is G = [cT^/(i + d)] ig , where is the pump 

pulse length and g is the gain per cm of the lasing medium. The gain 
* * 

g = KN n Q , where K is a constant, N is the density of excited alkali 
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atoms, and n Q is the density of the rare gas. In order to achieve maxi- 
mum gain, one would like the dead space d to be minimum. Thus our initial 
cell design had the windows of the cell in direct contact with the Li vapor, 
which will be at a temperature of about 1100°K. Conventional methods of 
sealing the windows, such as 0-rings could not be used. We investigated 
a metalization brazing technique, but although this could probably with- 
stand the temperature, it could not withstand the pressures. In addition, 

Li is extremely reactive at these high temperatures and very likely attacks 
both quartz and sapphire rapidly. 

In view of these difficulties, we analyzed the possibilities of using 
a heat pipe-type design in which the windows are cold, and isolated from 
the vapor by a zone of pure rare gas. It is then possible to use quartz 
windows with O-ring seals. Such a buffer zone, however, constitutes dead 
space in the optical cavity which reduces net gain. The minimum buffer 
zone is a function of the cell temperature, pressure, and diameter. If 
the temperature gradient between the center of the cell and the windows is 
too steep, Li will condense as suspended droplets rather than on the walls 
which will scatter the laser beam. Emp eric ally, the intensity of this 
’'raining*' appears to be directly proportional to the pressure, the temper- 
ature gradient in the buffer region, and inversely proportional to the 
cell inside diameter. From our experience with other heat pipes it seems 
that the temperature gradient should be limited to about T5°c/ (diameter) • 

(atmosphere) to avoid raining. The minimum cell diameter and total length 

♦ ^ 

are restricted by optical diffraction considerations; the minimum temperature 
gradient is limited by the allowable dead space in the cavity; and the gain 
per pass is determined by both pressure and total cell length. A simple 
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mathematical model incorporating all these effects was used to design the 
best practical cell geometry. The cell we have constructed will have an 
active zone of 20 cm, with a 10 cm buffer zone on each end, and an inner 
diameter of 0.6 cm. It consists of a stainless steel outer tube with 
flanges on each end. Quartz brews ter windows are sealed to the flanges 
with 0-rings. A slotted stainless steel tube with a stainless steel wick 
wrapped around it is inserted in the outer tube. The cell is in the pro- 
cess of being loaded and tested. 

To pump the Li-Xe we will use the second harmonic of the 1.338 |i line 
of a Nd:YAG laser. The Li resonance line is 6708 8 , but due to the Xe 
pressure broadening, the 6690 8 will be absorbed strongly. The calculated 
absorption depth is approximately 3 cm at 3 0$ Excitation of 3-9^ torr Li 
vapor pressure and the required pump power to achieve 30 $ excitation for 
100 ns in the above cell design is 3 kW. We will use a LilO^ doubling 
crystal and are presently assemblying the pump laser. 

Based on a report by York and Gallagher^ we have recently been examin- 
ing the properties of a high pressure Na-Na excimer. This system is dis- 
tinctly different from the low pressure Na^ molecular laser in which ground 
state Na^ molecules are directly pumped to excited states and lase to high 
lying ground vibrational levels. In the Na-Na excimer case ground state 
atoms are excited (as in alkali-rare gas excimers) and then combine with 

■K' 

ground state Na atoms to form Na^ (Na + Na -^Na^) in the presence of 
high pressure buffer gas. This process can be pumped much more efficiently 
because the ratio of Na/Na 2 is typically greater than 50 . In addition, the 
lower melting point of Na will reduce the raining and dead space problems. 
We are still investigating the Na-Na excimer as well as the other alkali- 
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alkali excimer systems. Our present cell design should be adaptable for use 
with any of the alkali-alkali systems. 
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C . Research Studies on Techniques for Frequency Conversion of High Power 
and High Energy Laser Radiation 
(D. M. Bloom, J. F. Young, and S. E. Harris) 

For the past period our efforts have been directed mainly at replac- 
ing the xenon formerly used as the phasematching gas with a more disper- 
sive vapor. Experiments with a two -metal heat pipe containing sodium and 
mercury indicate that mercury is not a suitable metal to be used in a two- 
metal heat pipe with the alkali metals. The failure of the sodium-mercury 
heat pipe is attributed to the alloying of the metals which greatly alters 
their vapor pressures. 

In order to understand, the failure of the sodium-mercury heat pipe it 

is necessary to understand the general operation of a two-metal heat pipe 

oven. Two-metal heat pipe ovens were first proposed and demonstrated by 

Hessel and Jankowski^ and later modified to allow arbitrary mixtures by 
2 

Vidal and Hessel. Detailed descriptions of the operation of this type of 
heat pipe oven are given in these papers. However, an important point not 
made in these papers and demonstrated in our experiments is that for proper 
operation the two fluids must be immiscible. Immiscibility is required to 
insure that the vapor pressures of the constituent fluids are not modified 
as they would be in a solution of the two. 

A study of binary alloys was carried out in order to find a suitable 
metal which was immiscible with sodium. Magnesium was found to be immiscible 
with sodium, strongly positively dispersive, and has a vapor pressure similar 
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enough to sodium to allow operation as a two-metal heat pipe with reasonable 
temperature gradients. A second two-metal heat pipe oven containing sodium 
and magnesium was then constructed. Experimental results with this cell 
indicate that this heat pipe oven is operating correctly and phasematching 
has been obtained to as many as 50 coherence lengths. This is the first 
phasematched system using two metal vapors to be demonstrated and the low 
phasematching ratio (N = ^ ) should allow operation at significantly 

higher sodium pressures. 

Currently, experiments are underway to determine the energy conversion 
efficiency obtainable with the present Na:Mg heat pipe oven and a modified 
version is being constructed in order to obtain a more homogeneous mixture. 

It is hoped that these systems will allow significantly higher conversion 
efficiencies than previously reported. 

Finally, construction of the birefringent filter is now complete and 
it will be installed in the laser oscillator at the conclusion of the ex- 
periments with the two-metal heat -pipes. Operation of the laser at 1.07^- Mm 
should allow higher conversion efficiencies in rubidium phasematched with 
xenon. 

REFERENCES 
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D. High Average Power Blue and UV Sources 

(G. W. Bekkers, J. F. Young, and S. E. Harris) 


In Appendix B of the July 197^- report it was demonstrated that in N a 
vapor a cw 9*^6 pm CO^ laser can be up-converted to 3321 & with a photon 
efficiency of 58$. The realization of this high conversion efficiency de- 
pends on tuning the pump frequency into exact two-photon resonance and also 
the incoming IR radiation has to approach resonance with the third level. 
This conversion efficiency is obtained in a single coherence length of 
metal vapor and is therefore quite broadband in both wavelength and aper- 
ture. However, the greatest advantage of this process lies in the possi- 
bility of high average power conversion. Our proposal is to use a 1.0790 pm 
Nd:YAL0^ laser which is exactly resonant with the 6s -7s transition in Cs and 
up-convert a 3 pn chemical laser line to U 556 S. 

The maximum power density which can be used in two -photon resonantly 
pimped processes is limited by two-photon absorption under the condition 
that W v J x = \/2 , where t is either the incident laser pulse length At 
or the decay time T^ of the second level, whichever is shorter. In Appen- 
dix A of the January 197^ report the last coherence length power conversion 
efficiency £ was calculated for the case that the coherence length is set 
by the upper level and the power density is limited by two -photon absorption 


e = 


2 ^ 

T&D 2 ^05 
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where the p’s are the dipole moments between levels indicated by the sub- 
scripts and tog is the half-power linewidth of the two-photon transition. 

£ is the power conversion efficiency of the incoming IR radiation into the 
generated UV or blue radiation. It is generally assumed that forbidden 
transitions like s-s or s-d are only lifetime broadened but dephasing col- 
lisions might cause pressure broadening. We were not able to find any 
literature on the broadening mechanism of these unallowed transitions. 

Therefore in our calculations we have to introduce an uncertainty range 

-b -1 

of = 5.3 X 10 cm in case of lifetime broadening to a maximum of 

Sci^ = 0.1 cm ^ in case of pressure broadening. Even in the latter case a 

power up-conversion efficiency £ = 200$ can be achieved. In case of third 

harmonic generation we can use the expression for £ given in the July 1974 

( 2 ) / 

report together with the condition J 1 ~ 1/2 . This leads to an ex- 

pected power conversion efficiency o f o.4$ in one coherence length in case 
of lifetime broadening. In case of maximum pressure broadening this num- 
ber is reduced by the square root -of the ratios of ■ . The result is 

a conversion efficiency of 0 . 029 $. 

A heat pipe is a device in which a chemical substance is evaporated in 
the center of a pipe and then condenses at the ends giving off its latent 
heat. The liquid substance is returned to the center by means of a wick* 

By filling the ends of the pipe with an inert gas the vapor pressure of 
the substance can be controlled. The great advantages of this device are 
the constant temperature along the pipe and the fact that the substance 
does not touch the windows. A stainless steel heat pipe using Cs vapor 
with argon as an inert gas is presently in operation at pressures from 
1-60 ram Hg. The wick is made from fine stainless steel mesh. 
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When crystals are inserted in a laser cavity the maximum average power 
of the doubled frequency can reach 2W. An example is a LilO^ crystal inserted 
into a Nd:YAG laser cavity running on the 0.9 ^6 line which delivers 0.5 W 
average power (Chromatix) . These systems can achieve a 40-50$ doubling ef- 
ficiency. Up-conversion in heat pipes can also achieve this conversion ef- 
ficiency and at the same time much higher average powers can be obtained. 

This is due to the high speed with which the alkaline atoms travel through the 
heat pipe. When atoms leave the evaporation zone, which is centered under 
the heating element, they enter the condensation zone. The nonlinear pro- 
cess only takes place in the evaporation zone where the density of atoms is 
constant. It has been determined that at the interface between evaporation 
zone and condensation zone the velocity of the atoms can be subsonic. For 
10 torr of Cs vapor the sonic speed is around 2^0 m/sec, This means that 
for a 1.5 cm long zone all atoms are replenished inbetween laser pulses up 
to high pulse repetition rates. Therefore if a high conversion efficiency 
can be achieved in one pulse, then high average power conversion is also 
possible, e.g,, a 1.5 cm long zone, a 0.5 cm laser beam diameter, a typical 
power density of 5 X 10 w/cm and a 50 nsec pulse leads to a IT W average 
power ^556 & radiation in Cs vapor at 1000 pps . The pulse length can even 
be much longer. 

At the moment the biggest unknown is the linewidth of the two- 

photon transition. This linewidth can be calculated by measuring two-photon 
absorption at several operating points, since all other quantities are ac- 
curately known. Once is known we can verify its value by measuring 

third harmonic generation. The next step is to up-convert a tunable 5 pm 
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source and verify the results with theory. When the outcome of this expert^ 
ment is successful then high average power conversion using a } (jm chemical 
laser should be attempted. 

We started out, not with a heat pipe, but with a 0.5 cm long stainless 
steel cell with sapphire windows at the ends and loaded with Cs . The cell 
was heated several times to 450°C without any side effects. However, when 
the Nd;YALO^ laser was propagated through the cell, the sapphire windows 
showed a black deposition around the location where the laser beam passed. 

The area of this deposition was 10-100 times larger than the beam area at 
that location. A telephone conversation with J. Reintjes at the "Naval Re- 
search Laboratory in Washington, who is measuring self-defocus ing in Cs 
vapor at 1.06 |_an, revealed that he also observed this deposition. These 
two independent observations make it unlikely that any contamination is 
involved. The physical origin of this deposition has not yet been pin- 
pointed. 

The discovery of the near coincidence of two Nd:YAL0^ photons at 
1.0795 an -d the 6s -7s transition in Cs led to a search for exact reso- 
nance. A smaller detuning from the second level means that a certain con- 
version efficiency can be achieved at a lower power density. Measuring the 
Nd:YAL0^ fluorescence spectrum showed that the desired 1.0790 wavelength 
was at the 4o $ level. By inserting an etalon into the laser cavity the 
Nd:YAL0^ laser now can be tuned over the range 1.0799 ->1.0789 (Jm, deliver- 
ing up to 5 per 50 nsec pulse. 

The first experiment was aimed at carefully measuring two-photon absorp- 
tion and calculating from these data. Two-photon absorption can be 
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distinguished from linear atomic absorption by its power dependence and from 
linear molecular absorption by its temperature dependence. At an atomic 

17 

density of 2 x 10 atoms per cc we measured a molecular Cs^ absorption of 
50 $ in a 1.5' cm long zone. When the power density was increased a new 
phenomenon occurred. The number of atoms in the volume of the laser beam 
is known and since only one-half of these atoms can be excited to the 7s 
level the maximum two -photon absorption is determined. However , the total 

v 

absorbed energy was about 50 times higher than this maximum amount and 
showed two peaks about 1.7 cm ^ wide and separated by 3.5 cm"'*'. From the 
observation of a threshold and the amount of loss we concluded this to be 
the result of a four-wave parametric oscillation; but a careful scan of 
the output spectrum did not reveal any of the two additional wavelengths. 
These two losses overshadowed the two-photon absorption, so that this ex- 
periment was inconclusive for . As a consequence we have to operate 

the Cs heat pipe at a few torr where the ratio of molecules to atoms is 
small and at power densities below the threshold of above described effect. 

While the Nd:YAL0^ laser was tuned over its 1.0799 -> 1.0789 pm range 
the third harmonic generation peaked at the expected 1.0790 pm wavelength. 
The measured value was 10 - 10 percent at the peak and the full width 

at half power was 1 cm \ This is wider than one would expect from a 0.3 
cm wide laser line. 

The experiment which is set up at the moment is to add the tunable 


output around 5 (jm of an optical parametric oscillator to the NdrYALO-, in- 

3 

put. The oscillator contains a *4- cm LiNbO^ crystal which is pumped by a 


doubled 1.0790 7m. This doubling is done in a CDA crystal. Tuning of the 
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oscillator is performed by changing the crystal temperature. The two beams 
are fed col linearly into the Cs heat pipe. Measurements of the up-conversion 
efficiency vs. wavelength are currently underway. 
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EFFECTS OF FOCUSING ON THIRD -ORDER NONLINEAR 

* 

PROCESSES IN ISOTROPIC MEDIA 

by 

Gary C. Bjorklund**' 

Microwave Laboratory 
Stanford University 
Stanford, California 


ABSTRACT 


A theoretical analysis of third-order nonlinear interactions of fo- 
cused laser beams is performed for the processes ^ ©^ , 

+ ©£ - ©^ -> ©^ , and ©^ - ©^ - ©,. -* ©^ • The total power and far- 
fleld beam profile of the generated radiation is related to the total 
powers of the fundamental beams, to the tightness and location of the 
focus, and to the value of the difference between the wave vectors of 
the generated radiation and driving polarization. The optimum degree 
of wave vector mismatch as a function of tightness and location of focus 
is determined for each of the three processes. The process ©^+©^-©^-» ©^ 
is found to be unique in that it is always optimized by focusing as tightly 
as possible* Experimental results, which verify the theory for the pro- 
cesses ©^ + ©g + ©^ and 00 ^ ^ ^ 9 are P res eflted. 


Work supported by the National Aeronautics and Space Administration. 
^Present address: Bell Telephone Laboratories, Holmdel, New Jersey 

07733 , (£01) 919-11917 . - 
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EFFECTS OF FOCUSING ON THIRD-ORDER NONLINEAR 


PROCESSES IN ISOTROPIC MEDIA 
by 

Gary C. Bjorklund 

Microwave Laboratory 
Stanford University 
Stanford, California 9^305 


I. INTRODUCTION 

Third-order nonlinear processes in isotropic media have recently been 
utilized for efficient frequency tripling of high power laser radiation , ^ 

Z )j r ^ 

for the production of tunable * and fixed frequency' 5 * coherent vacuum 

ultraviolet radiation, for up-conversion of infrared radiation, ^ for the 

production of tunable infrared radiation, ^ and for coherent anti-Stokes 

10-13 ' 

Raman spectroscopy. These processes involve frequency sum and/or 

difference mixing between as many as three different fundamental frequencies 
In order to achieve the high power densities required by these processes 
without damaging the windows of the cell containing the nonlinear medium, 
it is often necessary to focus the fundamental beams tightly enough to 
produce a depth of focus small in comparison to the cell length. Focusing 
in this manner dramatically alters the form of the functional dependence 
of the efficiencies on the mismatch between the wave vectors of the gener- 
ated radiation and driving polarization. 



This paper presents a theoretical and experimental investigation of the 
effects of focusing on the processes cd^ + + (u -> co^ , + <x > ^ -» 03^ , 

and -* cd^ . The optimum values of the wave vector mismatch, 

£k (defined as the wave vector of the generated radiation minus the wave 
vector of the driving polarization), are found to depend on the tightness and 
location of the focus. In the case of very tight focusing with the fundamen- 
tal beam waist region entirely contained within the cell, efficient generation 
can be achieved only for restricted values of . It is shown that 

must be non-zero and negative for the process -* and that 

must be non-zero and positive for the process -* 

Given a nonlinear medium with dispersion characteristics proper to produce 
the correct sign for Z& , each of these processes can be optimized by 
varying the tightness of focus. The process - eo -* is found 

to be unique in that it can always be optimized by varying the tightness of 
focus, regardless of the dispersion characteristics of the nonlinear medium. 

These focusing effects arise from the well known phase shift of 7T 
radians which occurs as a beam traverses the region of the focus. This 
phase shift also affects the driving polarization and, except for the pro- 
cess -* 0)^ , introduces a slip in phase between the driving 

polarization and the generated radiation. Destructive interference between 
portions of the generated radiation arising from different locations in the 
cell will result unless this slip in phase is compensated by a wave vector 
mismatch of the proper sign. 

In the following sections, a rigorous theoretical analysis is used to 
derive expressions which relate the total power and far-field beam profile 
of the generated radiation to the total powers of the fundamental beams, to 
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the tightness and location of focus of these beams relative to the dimen- 
sions of the cell, and to the value of the wave vector mismatch. Numerical 
evaluations of these expressions are presented for a wide variety of cases. 
Procedures for the optimization of the total generated power are discussed. 
Experimental results, which verify some aspects of the theory, are presented. 

It should be noted that the effects of focusing in the case of third 
harmonic generation in isotropic media have been previously considered by 
Miles and Harris and by Ward and. Hew. ^ The rigorous theoretical analysis 
of this paper is an extension of the work of Ward and New and of Kleinman, 

15 


et al. 



IX. THEORY 


la this section, expressions are developed for the total output power 
and beam shape of the generated radiation produced by the four -wave mixing 
processes + o)^ -> + cd^ - (x>^ -» , and - a> 2 - co^ -* 

The basic approach is to specify the spatial variation of the driving polar 
ization, perform a Fourier decomposition of the driving polarization into 
plane wave components, calculate the generated radiation field arising from 
each such component, and then to determine the total generated field by sum 
ming the contributions of all of these components. 

* The fundamental beams at ^ , and o> ‘are assumed to be lowest 

order Gaussian modes which propagate concentrically along the z axis with 
identical waist locations and identical confocal beam parameters. If the 
wave vectors of these beams in the nonlinear medium are taken as k, , k« , 

i c: 

and k^ respectively, and if the electric field amplitudes are taken as 
*10 ’ *20 , and respectively, then the fundamental electric field is 

given by 


i( 


ryt) = Re 


-i(o t -ico t -i(o t 

^(r) e + £„(?) e + £ (?) e ^ 


( 1 ) 


with , where n - 1,2, or 3 > given by the theory of Boyd and 



Gordon^ as 


^ n (*) = ^ n0 ex P (lk n z)(l + iS}' 1 exp 


- k n (x 2 + y 2 )/b(l + i£) 


. (2) 


In Eq. (2) ; b is the confocal parameter, defined by 


2t?t> 


b = 


0 


rw 2 

2noi 0 n 


2X 

0 , 2 

^ = 


(3) 


where is the beam waist radius, n is the index of refraction, 

is the vacuum wavelength, and Q is the far-field diffraction half angle. 
£ is a normalized coordinate along the z axis, defined as 


S = 


2(z -f) 


(*) 


where f is the position of the focus (beam waist) along the z axis. 

If the fundamental and generated beams are all assumed to be linearly 
polarized in the same direction, the vector nature of the electric fields 
as well as the tensor nature of the nonlinear susceptibility may be neg- 
lected, and the driving polarization of 03^ is given by 



) 


where 


M?) = 


| NX (- , co g , oi 5 ) Ej (r) Eg(r) E^r) 

(for + ODg + CD^ -» 0 ^) 

7 "X* 

| NX (- CDj^ ; m 1 , m 2 , -(Dj) E 1 {r) E 2 (r)E ? (r) 
(for ^ + m 2 - cd -»o^) 

| NX (- ^ ; oi 1 , -m £j -a> 3 ) E^r) E g (r) E^r) 
(for u> x - dig - ia -> ) 


(6) 


with X the nonlinear susceptibility per atom and N the number density 
of atoms of the nonlinear medium. [When two or three of the input fre- 
quencies are degenerate, the factor of ^ in Eq. (6) must be changed to 
| or J respectively.] 

The nonlinear medium is considered to occupy the semi-infinite space 
2 > 0 , with the space z < 0 considered to be a vacuum. Then from 
Eqs. (2) and (6), for the process cd^ + + <o , the driving polar- 
ization at is given by 


P 4 (7) . - HX (- <V <V ® 2 . <y E 10 E 20 E 30 • 


X exp (ik'z)(l + i t)~^ exp 


• " 2 + y 2 ) 1 


k'(x c 
b(l + 


it) 


B(«) 


(T) 


2.6 - 


where k 


f 


the wave vector of the driving polarization, is given by 


k' = + k 0 + = k 

\ d y 


and where 


B(z) = 


1 z > 0 

0 z < 0 


Similarly, for the process 


P^(?) = “NX (.a^;iD lf « 2f -® ) E^E q B(z) 


X exp (ik'z)(l + iP). 2 (l - i<?) ^ exp 


(- V." + iSk ' ) (x 2 + y 2 ' 

b (1 + F 2 ) 


where k' is now defined as 


k ~ k l + k 2 " k 3 9 


and k is defined as before* 


Finally, for the process , 


= - ■* (- «V«V ' V "V E 10 E 20 E 30 B(z) 


X exp (ik'z)(l + if) *(l-if) 2 exp 


(- k + ifk')(x 2 + y 2 ’ 
b(l + f 2 ) 


(8) 


(9) 


, ( 10 ) 


(ID 


, ( 12 ) 
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where k" is now defined as 


k ' = k l - k 2 - k 3 ’ 


( 13 ) 


and k is defined as before. 

The next step is to perform a Fourier decomposition of (t) to de- 
termine the amplitudes of the plane wave components of the driving polari- 
zation. , the amplitude of the plane wave component with wave vec- 
tor t , is defined by 




= (2tt) 


-3 


J * dx* J* &y" J* dz" P^ ) exp (- iif • /t”) 


(lh). 


, the generated radiation field arising from the plane wave 
component of the driving polarization with wave vector t , and P^j^y 
can be shown by Maxwell's Equations to be related by 

^X^xE^^fr) - k^E^Or) = exp (it . r) , (1?) 


where k^ and are the wave vectors of the generated radiation in the 

nonlinear medium and in vacuum respectively. Equation ( 15 ) may be solved 


subject to the boundary condition that the tangential component of 

vanish at the vacuum-medium interface. For K , K « K and 

x y z 

| - kj^ « | K*| , k^ , Kleinman^ has shown that 


(r) = i tSnk^z/k^] g (U^) exp (it * ?) P]^ 


( 16 ) 
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with 

%> = iz [K z - H + (K z + Ky)/2k 4 ] , (17) 

and 

1 

g(x) = [1 - exp (- x)]/x = f e~ xp dp . (18) 

0 

The wave vector mismatch, £k , is defined by the relation 

= k h - k' . ' ( 19 ) 

This choice of sign for £k is opposite from that of Ward and New 1 ^ but 
consistent with Miles and Harris.^" 

1 

The substitutions == k r + Zk. and p = - — b(£Z -£•)/ z may be per- 


formed in Eqs. (IT) and (l8) 

to yield (assuming Zk « , k' ) 


g(U K >) = ( 2 I) f dF ' eXp | 

lb (K -k'+ (K^ + K^)/ 2 k' - 2 k] (r -t)/z\ 

z x y i 

, (20) 

where £* - 2(z-f)/b and £ 

= 2f/b . 


The total generated electric field, E^(r) , is derived by combining 

the > 

CO 

03 00 


% (r) = J 

dK x / dK y / dK Z \(Z)& ’ 

(21) 

—CO 

-03 -00 


where E^(r) is defined by Eqs. (1 6 ), (20), and (lU). 
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r 


These integrals may be evaluated by first integrating over K z to 

yield the factor 6(z*-z) , then integrating over z* , and finally 

integrating over vi f y* j K , and K by means of the formula 

x y 


1 

y exp [- iqx - isx 2 ] dx = (tt/1s) 2 exp (iq 2 As) 


For £D^ + -> , the result is 


E^(r) » i - — 7Jk 2 bX (- E^E^E^ exp (ik'z)(l + if) 

2 H 


-x/2 2 ■» ~i o 

- k (x +y ) * 

b(l + i£) % 


X exp 


exp t- (ib/2) £k (e/ - ?)] 

(1 + tr ) 2 


while for -» , 


= i Tik^bX (- aV'V'V '“3) E 10 £ 20 E 30 exp (ik ' z) 


exp [- (ib/2) £k (P' - £)] f- (x 2 + y 2 ) 

df ' exp 

(l + ir)(k" - ik'f')H I bH 





and for - <o -» cd^ > 


3N 


\(?) = i *V X ( ' ^4 ;m i ' ~ (X> 2 > •° ) 3 ) E 10 E 20 E 30 eXp (ik ' z) 

2k 


exp [ - (ib/2) (£' - £)] 

• x I d£' exp 


f 


(1 - ±&'){k" -ik'8')H . 


, 2 2 S 
- (x + y ) 


bH 


(25) 


k' is defined by Eqs, (8), (11), and (13) in Eqs . (23), (24), and (25) 
respectively, k" is always given by Eq. (8) and the function H is 
defined by 


H b H (£,£') = 


a + g ,g ) 

(k" - ik'f') 


- 1 


(r - s) 

k / 


(26) 


If the nonlinear medium is contained in a cell of length L with in- 

, put window located at z = 0 and with the location of the focus at z = f , 

then Eqs. (23), (24), and (25), with £ = 2(L-f)/b , z = L , and t = 2f/b , 

completely describe the amplitude of the generated field at the plane of the 

output window. The functional dependence of the field distribution upon the 

transverse coordinates x and y scales with’the value of k' , b , and £ . 

It can be seen from Eq. (23) that for the process -» , 

E^C?) is always a lowest order Gaussian mode with the same confocal parameter 

and waist location as the input beams. However Eqs, (24) and (25) 

show that for the processes u) + oj - w -+u>. and u> - w - 

12 3^ 1 2 

W 3’" > " t0 4 * is n °t l n general a single Gaussian mode except 

when k" = k ? . The multimode output of the difference 
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mixing processes may be explained as resulting from a mismatch between the 
size of the driving polarization beam waist and the size of the generated 
radiation beam waist which would correspond to the confocal parameter of 
the input beams. This mismatch may be represented by the ratio k" fk' • 
When k'* /k' =1 , these beam waists are of the same size, and a lowest 
order Gaussian beam with the same confocal parameter and waist location as 
the input beams can be produced. The sum mixing process -» 

always has /k' = 1 and thus always produces a lowest order Gaussian 
beam output. 

For all three mixing processes it may be seen that E^(r*) is cir- 
cularly symmetrical about the z axis and that the total generated power 
at is thus obtained by performing integrals of the form 


J 2ttR |e(R*)| 2 dR , (27) 

• 0 

2 2 “ 

where R = (x + y ) 2 . It is convenient to define the dimensionless func- 

b f k^ 

tions Fj (b^k , — , — , j—r) by the relations 




b X l E xo E 20 E 30^ 


X f 27TR |e(R )| 2 dR , (28) 

where j = 1,2, or 3 implies ^ + U3 2 + <D^ 0^ , C0 2 + <U 2 - £D -» , or 

- (D . Although the spatial distribution of E^(r*) depends 

upon the value and units of k' and b' , the integrated intensity repre- 
sented by the Fj depends only on the dimensionless parameters b/Sk , b/L , 
f/h , and k" /k' . 
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-1 


With all wave vectors expressed in units of cm , b in cm, N in 

, 3 

atoms/cm , X in e,s,u t per atom, and with total power in watts, the 
total generated power, P^ , is given by 


->k k i k - 


= (6.318 X 10'") ? • F j , - , - , — j , (29) 


where P^ , Pg , and P^ are the total powers of each of the fundamental 


beams . 


For the process cd^ + 03^ + ^ , 


/ b f k" \ 

^ bZk , - , - , — 

\ L L k' / 


f d£ 

4 


exp [- (ib/2) £'] 

(1 + if /) 2 


• (30) 


Thus F^ (b £k , b/L , f/L , 1 ) is equivalent to the quantity |l(3,£k,£,£)|‘ 

1^ 

as defined by Ward and New (who define with the opposite sign). In 

the tight focusing limit where # b is short compared to the cell, and the 
entire focal region is contained within the cell (£ , £-*<») } the integral 
may be solved in closed form to yield 


7T 2 (b£k) 2 exp [b£k/2) for £k < 0 


^^, 0 , 0 . 5 ,!) - 


0 


. (31) 


for £k ^ 0 




2rrRdR 


with 

( 11 ) 


For the process + cd^ - o>^ , 

j b f k" \ 2k" 

Fp ( ^ ~ * j I = / 

2 \ L L k" / 71b J n 


^ exp [** (ib/2) £k S'] 
i d£' 

4 


(l + iS')(k // -ik'S')H 


X exp 


- R 


bH 


(32) 


H defined by Eq. (26) and k' and Vl" defined by Eqs. (8) and 
In the tight focusing limit and in the case k^ = k ' , 


F 2 (bAk ,0,0.5,!) = 7 T exp [- b\£k\ ] 


(33) 


For the process - o> — , 

2k' 

7lb 


/ b f k tf \ 2k' Z. 

{ b£k , - , - , — ) = / 

\ L L k' / 


2rrRdR 


Z. exp [- (ib/2) S'] 

/ dS' 

(1 - iS' 


(1 -iS')(k" -ik'g')H 


X exp 


- R 


bH 


2 n 


( 5*0 
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with H defined by Eq. (26) and k' and k* defined by Eqs * (8) and 
(1,5). In the tight focusing limit and in the case k* = k' , 


0 


for £ o 


F 3 (b£k, 0,0.5, 1) 


7 T 2 (b^) 2 exp [b^k/2] for £k > 0 



III. EVALUATION OF THE 


The functions were numerically evaluated as continuous functions 

of b/ik for various sets of discrete values of b/L } f/L , and k Jk' . 

From Eqs. ($0), (32)j and (3*0 it is evident that the evaluation of each 

F^ involves an integration over £■' which may be put in the form of a 

Fourier transformation. An efficient discrete Fast Fourier Transform rou- 
19 

tine 7 was employed to numerically evaluate the F^ on an IBM computer. 

This routine is capable of performing a ^09^ element complex discrete 
Fourier Transform in less than *1- seconds of computer time. One such 
transformation is sufficient to completely specify as a function of 

for one set of values of b/L and f/L , or to evaluate |Ej^(if)|^ 
as a function of b^k for one set of values of b/L , f/L , k fk" , and 
R . Fg and were calculated as functions of b£k for one set of 

values of b/L , f/L , and k" /k' by evaluating |e^(r)|^ for *1-0 discrete 
values of H and then integrating over R by means of the trapezoidal rule 
Thus the determination of Fg or F^ as a function of b£k for a single 
set of values of b/L , f/L , and k" /k' required the Fast Fourier Trans- 
form routine to be run *1-0 times and required approximately 3 minutes of 
computer time. 

. . \ 

The behavior of the F^ near the tight focusing limit was studied in 
detail. It was determined that, when b/L ^ 0.1 and when b ^ f £ L-b , 
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the F j as functions of bZsk may, to a goo 4 approximation, be represented 
by the curve for F vs. b£k corresponding to the infinitely tight focus- 
ing limit where b/L =0.0 and f/L =0.5 • For instance, the maximum 

deviation of the numerically computed F^ (bZk , 0 . 1 , 0 . 5 , l) from the 
analytically derived F (b£k , 0 , 0.5 , 1) of Eq. (jl) is less than 0.1 
for all values of bZk . 

F- for b/L £ 0.1 and f/L = 0.5 is plotted in Fig, 1; F 0 for 
b/h s 0,1 , f/L = 0.5 , and k /k' = 1 , 1.1 , 1.5 , 2.0 , and 3,0 are plotted 
in Fig. 2; and F^ for b/L £ 0.1 , f/L = 0.5 , and k"/k' = 1 , 1.1 , 1.5 , 
2,0 and 3*0 are plotted in Fig. 3 . The analytical results of Eqs, (31); 

(55 )> and ( 35 ) were used to supply the curves for b/L =0 , while the 

numerical results were used for b/L / 0 . 

In the case of tight focusing, the F_. , when viewed as functions of 

b£k , have the mutual characteristic of possessing a single large peak. As 
the value of b£k departs from the optimum value, the F_. approach zero in 

a nearly asymptotic manner. In the case of b/L slightly more than zero, 

some sidelobe structure appears far in the wings. However, as long as the 
fundamental beam waist region is contained entirely in the cell, these side- 
lobes are always less than 0.1$ of the magnitude of the peak value of F 

j 

(This is in distinct contrast to the plane wave case, where sidelobes to the 
main peak occur with magnitudes as high as 0 f the maximum value. ) The 
magnitudes and locations of the peaks of the F^ are given in Table I. 

In the case of very loose focusing, the F^ approach the plane wave 

limit, as given by Miles and Harris,^ of 

/ b k" \ hi? /ZkL 

lim F b£k , ~ , 0.5 , — I = —5- sine ( 

b/L -> ~ 3 \ L k' / b \ 2 
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TABLE I 



b£k for tight focusing cases 


F. 

J 

Peak 

Value 

Value of 
bZik at 
peak 

: Range of b£k 
1 for which 
F. ^ 50$ of peak 

Range of bZ)k 

for which 

F. ^ 10 $ of peak 

Range of b£k 

for which 

F, ^ 0.1$ of peak 







? l ( b ^ , 

0 , 0 , 

•5 , 1) 

50 

- 2,0 

- 

4.1 


- 

0.8 

- 

6.7 

to 

- 0.3 

- 

12.5 

to 

- 

0.1 

F x (b£k, 

0.1 , 

0.5, 1) 

5.4 

- 2.1 

- 

4.1 

to 

- 

0.8 

- 

6.7 

to 

- 0.3 

- 

12.5 

to 

+ 

0.1 

F x (b£k, 

0.1 , 

0.75, 1) 

5.5 

- 2.0 

1 

4.1 

to 

- 

0.9 

- 

6.7 

to 

- 0.3 

- 

12.5 

to 

+ 

o .6 

F 1 (b^c , 

0.1 , 

0.85, 1) 

5*4 

- 2.2 

i ~ 

4.1 

to 

- 

0.9 

- 

6.7 

to 

- 0.3 

- 

12.4 

to 

+ 

1.9 

F 2 (bAk, 

0 , 0 . 

.5, 1) 

9.9 

0.0 

i 

0.7 

to 

+ 

0.7 

- 

2.3 

to 

+ 2.3 

1 - 

6.9 

to 

+ 

6.9 

F 2 (bac, 

0.1 , 

0.5 , 1) 

8.5 

0.0 

• 

li 

0.8 

to 

+ 

0.8 

- 

2.4 

to 

+ 2.4 

, - 

7.0 

to 

+ 

7.0 

F 2 (b^k, 

0.1 , 

0.5, 1.1) 

7.7 

0.0 

I - 

0.8 

to 

+ 

0.8 

- 

2.4 

to 

+ 2.4 

j ^ 

7.0 

to 

+ 

7.0 

F 2 (b^, 

0.1 , 

0.5, 1.5) 

5.* 

0.0 

- 

0.8 

to 

+ 

0.9 

- 

2.4 

to 

+ 2.7 

l 

7.0 

to 

+ 

8.7 

F 2 (b£k, 

0.1 , 

0.5 , 2 ) 

3.8 

0.0 

- 

0.8 

to 

+ 

0.9 

- 

2.b 

to 

+•3.3 

- 

7.0 

to 

+ 

11.4 

F 2 (bA, 

0.1 , 

0.5 , 3) 

2.1 

0.0 

- 

0.8 

to 

+ 

1.0 

- 

2.4 

to 

+ 4.6 

- 

7.0 

to 

+ 

16.9 

(b£k , 

0 , 0 , 

•5,1) I 

5*3 

1 

+ 

0.8 

to 

+ 

4.1 

+ 

0.3 

to 

+ 6.7 

+ 

0.1 

to 

+ 

12.5 

F ? (b£k, 

0.1 , 

0 . 5 , 1 ) 

5.3 

1 

+ 

0.8 

to 

+ 

4.1 

+ 

0.3 

to 

+ 6.7 i 

+ 

0.1 

to 

+ 

12.6 

F 5 (b<3c, 

0.1 , 

0 . 5 , l.i) ! 

4.8 

+ 2.1 

+ 

0.8 

to 

+ 


+ 

0.3 

to 

+ 6.7 ! 

+ 

0.1 

to 

+ 

12.6 

F 3 (b£k, 

0.1 , 

0.5 , 1.5) 

3.4 

+ 2.1 | 

+ 

0.8 

to 

+ 

4.2 

+ 

0.3 

to 

+ 7.1 

+ 

0.1 

to 

+ 

13.8 

(b£k , 

0.1 , 

0 . 5 , 2 ) 

2.3 

+ 2.1 

+ 

0.7 

to 

4- 

4.4 

+ 

0.3 

to 

+ 7.8 

+ 

0.1 

to 

+ 

16.3 

F 5 (b£k, 

0,1 , 

0.5 , 3) 

1.4 

+ 2.1 

+ 

0.7 

to 

+ 

4.9 

+ 

0.3 

to 

+ 9.3 

+ • 

0.1 

to 

+ 

22.0 













Families of curves of F^ vs, b£k for various locations of the 


The value of k /k is irrelevant to the behavior of the F. xn the plane 

J 

2 . 

wave case. The peak of this sine function is offset from zero by £k = 

for F^ and by £k - h/b for F^ , while for F^ the peak is always at 

zero. 

The complete evolution of the form of the from the tight focusing 

limit to the plane wave limit is demonstrated in Figs, 4 and 5 for the case 
of F x 

focus (various values of f/L ) are plotted in Figs. 6, 7j 8* 9> an£ * 10, 

where each figure corresponds to a different, constant value of b/h 

Figure 6 illustrates that, for tight focusing, vs. b£k is relatively 

unchanged as long as the fundamental beam waist region is entirely contained 

within the cell. Figure 10 illustrates that, for b much larger than L , 

F. • vs . b£k is relatively unchanged as long as the cell is located close 

to the beam waist location. Figures 7j 8, and 9 show the more complicated 

behavior of the intermediate cases. In every case, however, when the beam 

waist location is far outside of the cell, the location of the peak of the 

F. lies very near to b£k ~ 0 . 

3 

The dependence of the mode structure of the generated radiation upon 
bAk may be determined by monitoring (E^(R)j as- a function of the trans- 
verse coordinate R for each value of b£k . For the processes “l +CD 2-“' 
and ^ , when k" /k' ^ 1 , the far-field intensity distri- 
bution evolves from a nearly Gaussian shape for bAk ^ 0 , to a circularly 

symmetric ring pattern as b£k becomes large and positive. This nearly 
:Gaussian shape is in close agreement with the far-field intensity distri- 
bution which would arise from a lowest order Gaussian mode with the same 
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0.5 and b/L =0,1,2, and 5 . 
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Fig. 6— F 7 vs. bAk for b /l = 0.1 and f /h = 0.5 , 0.75 > 0.85 , 1.0 , and 1.5 . 
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Fig. 10 — vs. b<ik for b/L = 10 and f /l ~0.5 , 1.0 , 1.'5 j and 2.0 . 


wavelength as the generated radiation and with the same confoeal parameter 

and waist location as the fundamental beams. The circularly symmetric 

19 

ring pattern is of the type which has been observed by Terhune and others 
for anti-Stokes radiation produced by three-wave mixing. The results of 
this section thus indicate that this ring pattern will appear only when 
b£k is positive. A typical evolution of the far-field pattern as £k 
is changed is presented in Fig. 11 for the process + ^2 ~ ^ 
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Fig, 11 — Far-field intensity distribution for several values of bAk for the process 
© 1 + ~ k* flu' == 5.0 • The distribution in each case is' 

calculated at the position of the output window of a cell with b/h = 0.1 
and f/L =! 0.5 * 



IV* A PHYSICAL BASIS FOR THE BEHAVIOR OF THE 'F 

Several aspects of the behavior of the F. , including the optimum 

J 

values of b£k and the invariance with respect to the exact location of 
the focus in the tight focusing case, may be explained by a simple model 
which takes into account the additional phase shift caused by focusing. 

From Eqs , (2) and (6), it can be seen that the amplitude of the driving 
polarization strongly peaks in the region of the fundamental beam waists. 

A region of significant generation may somewhat arbitrarily be defined as 
bounded by the z values at which the driving polarization falls to 10# 
of its peak value. In the case of tight focusing, when b is much less 
than L and the fundamental beam waist region is located entirely within 
the cell, the boundaries of the generating region are located at £ - - 2 
and £ = .+ 2 (or z = f-b and z ~ f + b ). In the case of loose focus- 
ing, when b is much larger than L , the boundaries are located at the 
input and output windows of the cell. 

.From Eq. (2), it can be seen from the term 

i 

(1 + i£) ^ = (1 + & 2 ) 2 exp (- i 'tan * 6) 

that a focused lowest order Gaussian beam undergoes a shift in phase given 
by tan"^ £ as it propagates through the waist region. The driving 
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polarization given by Eq. (6} then experiences a phase shift of 3 ban 1 £ 
for ^ + ci^ -»a^ , tan" 1 £ for o> 1 + - a> , and - tan" 1 £ 

for - 0 )^ - 03^ . If it is assumed that the optimum conversion ef- 

ficiency occurs when the generated radiation is nearly a lowest order Gaussian 
mode with the same confocal parameter as the fundamental beams, then the gen- 
erated radiation may be considered to experience a phase shift of tan 1 £ 

Thus a slip in phase between the driving polarization and the generated rad- 
iation of 2 tan" 1 £ for 0^+0^ + cd^ -> o ^ and of - 2 tan" 1 £ for 
^ ^ ^ occurs, while no slip in phase occurs for - co^ -» oo^ . 

This slip in phase causes destructive interference between radiation gener- 
ated in different portions of the generating region. * 

The slip in phase due to focusing may be compensated to some degree by 
the wave vector mismatch, , introduced by the dispersion of the nonlinear 

medium. The optimum value of £k , ^ opt * is that value which exactly 
cancels the total slip in phase which occurs as the generating region is 
traversed. 

In the case of tight focusing, when the entire generating region is con- 
tained within the cell, this requirement- yields 




opt 


for 

“i 

+ 

“2 

+ 0 ), 
3 



for 


+ 

°’ 2 

- m 3 


f 

for 

ox 


03 ^ 


-> 0 ), 



1 


2 

3 

4 



( 37 ) 


in good agreement with the exact results. Since significant generation oc- 
curs only in the generating region, it follows that tfcis behavior should be 
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independent of the exact location of the focus so long as the entire gen- 
erating region is contained within the cell. 

In the case of loose focusing, when the entire cell is located in the 
beam waist, these arguments yield 



®1 + ^ 


( 58 ) 


again in good agreement with the exact results . 


I 


V. OPTIMIZATION OF THE TIGHT FOCUSING CASE 


In this section, procedures for the optimization of the total gen- 
erated power are developed for the case of tight focusing (b/L ^ 0.1) 

with the fundamental beam waist region entirely contained within the cell. 

1 

If limiting processes such as' pump depletion, absorption, breakdown, satur- 
ation, thermal defocusing, and the quadratic Kerr effect are neglected, then 
the total generated power at is always. given by Eq. (^1) and thus obeys 

the proportionality relation 

P !t “ P 1 P 2 P 3 n2>c2 Fj(b4c,b/L, f/L,k"/k') . (39) 


>idered to be constant 


auu 


The parameters ^ ’ anc * fa' are cons 

thus the optimization procedure reduces to maximizing the quantity 

(b^k , b/L , f/L , k" fa' ) by varying the parameters N , b£k , b/L , 
and f/L 


When N is a parameter which is independent of b£k , b/L , and f/L , 
the optimization procedure further simplifies to increasing N to the • 
highest possible value while independently maximizing F . By the re- 
sults of Section III, the maximization conditions for the F. under these 
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given by 


tight focusing conditions are then simply to set bAk to bAk 


opt 


( 2 

for 


+ U) 2 

+ 0 ), -> 03 . 

3 

bAk ^ < 0 

for 

03 - 

+ a>_ 

- m, 

opt J 


1 

2 

3 

1+2 

for 

“1 

" m 2 

- 03 ^ -> 03 , 

regardless of the exact values of 

b/L 

and 

f/L 

(so long 


(Uo) 


and b^f^L-b ) , 

The parameter bAk may be adjusted to its optimum value by means of 
altering b or by altering £k . The value of b may be straightfor- 
wardly altered by varying the tightness of focus, but in order for the tight 
focusing approximation of this section to remain valid, b/b must remain 
less than 0.1 ♦ In order to make Ak a parameter which is independent of 
N , another medium with negligible nonlinearity but appreciable dispersion 
may be mixed with the nonlinear medium. ^ Ak could then be controlled by 
varying the concentration of the added medium. The process + a) o - <13^ *-> 0)^ 
is unique in that the optimum value of b is always zero, regardless of the 
value of Ak . It should be noted that the procedure of optimizing by ad- 
justing b is valid only if the nonlinear medium can tolerate the high power 
densities introduced by very tight focusing. If this is not possible, the 
best compromise is to reduce b to the smallest permissible value and then 

to adjust Ak to reach bAk , 

opt 

The value of F_. is severely reduced if b£k deviates from its optimum 

value. The tolerances on the value of bAk which correspond to F. main- 

J 

taining at least 10$, and 0.1$ of its peak value are given in Table I 
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of Section III. It is particularly important to note that the process 
03^ + tog + -»a^ requires that b£k be non-zero and positive while the 

process - cn^ - -> 0 )^ requires that b^k be non- zero and negative. 

The generated total power at cu^ under these optimized conditions is given 
by Eq. (29) with F. set to its peak value from Table I. 

Different optimization procedures are required when N ig a parameter 
which is not independent of b£k , One particularly important case occurs 
when £k is constrained to be proportional to N while b is constrained 
to be constant and N is a free parameter* In this case, the quantity which 
must be maximized in the optimization procedure is £k (b£k , b/L , f/L , 
k." fk.' ) . This quantity is not dimensionless, and depends on the absolute 

value of b Therefore it is appropriate to define the dimensionless 
quantity . - 


G,. (b£k , b/L , f/L , k"/k' ) = (£kb) 2 F^ (b£k , b/L , f/L , k" /k ' ) , (bl) 


and to realize that the quantity which must be maximized is then 

(l/b ) 2 Gj (b£k , b/L , f/L , k"/k' ) . 

is plotted for b/L = 0 and b/L = 1 with f/L = 0*5 > 0.75 , 

and 0.85 in Fig. 12. The location and magnitude of the peaks of G ^ are 

summarized in Table II. - Here it is apparent that, as in the case of the 

F, , the value of b/L in the G. may be to a good approximation be re- 
— 3 

placed by zero for b/L less than 0.1. In addition, for tight focusing, 
the Gj are independent of the value of f/L as long as the entire input 
beam waist region is contained within the cell. G~, is plotted far f/L = 0.5 , 
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. 1 , 0 . 5 ,!) 

G^bAk 

7 0 
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7 0 
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Gg(b£k 

70 

.1 , 0.5 , 3 *o) 
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7 0 
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TABLE II 


'/L, f/L, k7k') 


for tight focusing cases 




Peak Value (s ) 

Value of bAk at Peak 



46,3 

- 4.0 

46.3 

- 4.o 

^5.9 

- 5.8 

V7.2 

- 4.1 

!)6.8 

- 5.7 

53. ** 

- 2.0 

53 + 

+ 2.0 

53 + 

- 1.8 

53- 11 

4 - 1.8 

48.5 

- 1.8 

49.0 

+ 1.8 

34.2 

- 1.8 

4o.6 

+ 2.6 

23.8 

- 1.8 

4o.3 

+ 3.0 

13.3 

- 1.8 

45.0 

+ 5.0 

46.3 

+ 4.o 

45.7 

3.8 

4l.5 

3.8 

30.8 ' 

4,4 

23.6 

\6 

1 6.6 

5*0 
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b/h = 0.1 in Figs. Ij5 and 1^. is plotted for f/L = 0.5 t b/L = 0.1 

in Fig. 15 . 

Since N is the only free parameter, the optimization procedure is 
to adjust b£k ■ to maximize by adjusting N . Let (b£k)^ t ^e t * ie 

value of b£k as given in Table II which corresponds to the peak of • 

Define the coefficient of proportionality, (X , by the relation = ON . 
Then the optimization procedure reduces to adjusting N to the value 
N . = (b^k) fab . 

opt v 'opt' 

From Eqs . ( 29 ) and (4l), the total generated power in this optimized 
case is 


( 6.318 x 10 " 4 ) 


k Q k l k 2 k 3 

k^k' 


X%*,P P G. 

1 2 3 1 


a 




(U2) 


where is given its peak value from Table II. 

It may be noted that in order to reach this optimized condition, a 
must be negative for the process + ^ ^ t while d must be 

positive for the process -> 03^ . The process -> 

is unique in that Gg possesses peaks for positive and negative and 

thus 0! may be either negative or positive. 
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Fig. 13 --Gg vs. b£k for b/L ^ 0.1 ; f/'L = 0.5 ; and k"/k' = 1.0 , 1.5 > and 3.0 




G\ 

\jJ 



Fig. 1 vs. b£k for b/b ^ 0.1 ; f/L = 0.5 » ant * k A* - 1*0 > 1*1 ? 1 *5 > 2.0 , and 3*0 



VI. experimental results 


An experiment has been conducted which verifies the theory of this 
paper for the process + 0)^ 3^ (third harmonic generation) in 

the case of very tight focusing near the center of the cell. In these ex- 
periments the nonlinear medium was Xe gas, the fundamental beam had a wave- 
length of 35 J ;7 and the generated radiation had a wavelength of 1182 
For this process, Xe possesses both a high nonlinear susceptibility and 

the proper dispersion characteristics to produce a negative value for Ak. 

6 

The experimental set-up was identical with that of Rung, et al,, ex- 
cept that the 35^T & beam was passed through a lens -diamond pinhole -shearing 
aperture spatial filter in order to produce a very high quality diffraction 
limited beam with well known parameters. 

In this experiment, approximately 0.1# of the 35*4-7 & radiation was 
converted to 1182 & radiation. The fundamental beam was focused at the 
center of the cell and b/L was 0.025. Figure l6 shows the relative power 
output at 1182 & vs. Xe pressure with no Ar in the cell. The agreement of 
the triangular data points to the theoretical function G^(bAk, 0, 0*5 , 1) 
is excellent. The departure of the circular data points from theory is 
caused by the multimode character of the fundamental beam which results 
from an incorrect shearing aperture diameter. The process of fitting 
to the data yields the coefficient of proportionality between Ak and the 
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RELATIVE POWER OUTPUT AT 1182 A 


THEORY- 



Xe PRESSURE (Torr) 


Lg. l6— Relative output at 1182 X, vs. he pressure* (The triangular points represent 
data from a run in which the shearing aperture was set to the proper diameter 
The circular points represent data from a run in which the shearing aperture 
was open tpc far. The solid curve represents the best fit to the function 
G, (b~k , 0 , 0*5 ) -) to the triangular points.) 
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atomic number density of the nonlinear medium, and thus provides an alter- . 
native method of determining the dispersion between the fundamental and 

-IT 

generated wavelengths. The fit in Fig. 1 6 yielded 2k = (- 5*99 X 10 ^Xe 
(cgs units). 

Some aspects of the theory of this paper for the processes 03^ +03^ 

and 03 +03, + 03 -»03, have been verified in a series of experiments per- 

1 1 3 4 • 

k 

formed. by Rung. These experiments concerned the mixing of a fixed frequency 

pump at 2660 & (03 = 03 ) with a variable frequency signal or idler 

P 

(cn ,03. = 03 ) in Xe gas to produce tunable coherent vacuum ultraviolet 
radiation. Rung 1 s experimental configuration involved very tight focusing 
of the fundamental beams. Figure 17 shows the maximum possible spectral 
range which would be predicted for each process on the basis of simple 
addition or subtraction of the available fundamental photon energies. 

Figure 17 also shows the range of tunable VUV output actually observed by 
Rung. The process 03^ + 03^ + 03^ -»03^ produced a detectable signal at 03^ 
only when was near to, and on the short wavelength side of a strong Xe 

line. This is a necessary condition to produce anomalous dispersion and thus 
the negative values for 2k necessary to optimize G^ . The process 
03 1 + 03 - 03 ->03i , however, produced a detectable signal over the entire 

maximum possible spectral range. This result would be expected since G^ 
may be optimized for either negative or positive 2k . Thus the theory 
qualitatively explains the results of these experiments. 


- 66 - 



TUNABLE VUV 


(Uj = 7200 A — ► 5320 A 
oj s * 5320 A — «► 4200 A 
(iip 3 2660A 



' ^ 

Fig, 17 --Summary of the tunable VUV output achieved by Kung, = 

pump frequency at 266 0 8 , = signal frequency tunable 

from 5320 8 to U200 8 , <0^ = idler frequency tunable from 

7200 8 to 5320 8 . . 





VII. CONCLUSIONS 


This paper has presented a theoretical and experimental investigation, 
of the effects of focusing on third-order nonlinear processes in isotropic 
media. A rigorous theory which relates the total power and far-field beam 
profile of the generated radiation to the various parameters of; the funda- 
mental beams and to the wave ‘vector mismatch has been developed. In the 
case of tight focusing with the fundamental beam waist region entirely con- 
tained within the cell, it has been shown that the process ->£0^ 

can be efficient only in media whose dispersion characteristics produce" a 
negative wave vector mismatch and that the process can 

be efficient only in media whose dispersion characteristics produce a posi- 
tive wave vector mismatch. The process - - in can be efficient 

. -■ -r ■ , . i . •„ • - , . . § # . . . 1 

in media with all types of dispersion characteristics and is always optimized 
by focusing as tightly as possible. 

Hie improved understanding contributed by this paper of the far-field 
beam profile of the generated radiation for the process ->cn^ 

has important .applications for coherent anti-Stokes Raman spectroscopy. By 
controlling the value of b£k , the mode of the output radiation may be 
varied from near Gaussian to a ring shaped pattern of variable cone angle. 

The cone angle of this ring pattern may be increased to the point where the 
generated radiation and the Gaussian fundamental beams become completely 



spatially separated in the far field. The advantages in detectivity which 
accrue from spatial separation between the generated radiation and funda- 
mental beams may thus be secured for geometries which involve colinear 
propagation of the fundamental beams. 

In the VUV and soft x-ray spectral regions index of refraction data is 


often completely lacking and even the sign of the wave vector mismatch pro- 
duced by various media may be unknown. In this case, if £k is of the 
wrong sign and/or if b is of the wrong value, the efficiencies of the 
processes and ^l “ ^2 " ^3 raay be re< ^ uce< ^ by 

several orders of magnitude. Thus the process (0- + 0) - cn ->cd. with very 

1 £ 3 ^ 

tight focusing has the highest probability of successfully producing a suf- 
ficient output for detection. 
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